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ABSTRACT 

Two-dimensional systems of C20 fullerenes connected to each other by strong covalent bonds 
have been investigated. Several isomers differing in the type of intercluster bonds have been 
revealed. The hfetimes r of the (C2o)mxM complexes with M = 2 and 3 at T = (1800 — 3300) 
have been directly calculated using the molecular dynamics method. It has been shown that 
these complexes lose their periodic cluster structure due usually to the coalescence of two 
or several neighboring C20 fullerenes. The activation energy of this process determined by 
analyzing the r (T) dependence appears to be Ea ~ 2.5 eV in agreement with the calculations of 
the heights of the potential barriers preventing the coalescence. At high temperatures T > 2400 
K, the decay of C20 fullerenes entering into the complex is possible. 



PACS: 36.40. Qv,61.48+c,71.15Pd 



The minimum possible fuUerene C20 has attracted considerable interest since its discovery 
in 2000 [1]. Great attention is focused on the study of the possibility of existing macroscopic 
crystaUine forms based on C20 fuUerenes similar to Cgo fuUerites [2,3]. According to theoretical 
works [4-7], C20 clusters in crystals are connected with each other by strong covalent bonds. 
This property distinguishes the C20 fuUerite from the Ceo fuUerite, where the bond between Ceo 
clusters is due to a weak van der Waals attraction. Interest in the C20 fuUerite is stimulated 
in particular by the expectation that the superconducting transition temperature in it will be 
higher than that in the Ceo fuUerite (significant enhancement of the electron-phonon interaction 
with a decrease in the size of fuUerites from Ceo to Cse, C28, and C20 was predicted 



There is experimental evidence for the existence of (C2o)f^ complexes with N — 3 13 
[9] (which are likely chains of C20 fuUerenes [10, 11]) and bulk phases [12,13] in which C20 
fuUerenes form a crystal lattice as atoms in normal solids. The structure of these phases is still 
unknown (the transmission electron microscopy data [12] are inconsistent with the electron 
diffraction data [13], indicating that the phase observed in [12] differs from that observed in 
[13]). The density functional theory investigations [4,5,7] show that it can be both simple 
cubic and body-centered (cubic, tetragonal, or orthorhombic) . For the same lattice type, the 
intercluster bonds can have various forms [7] , leading to a large variety of atomic configurations, 
candidates for the C20 fuUerite. All of these configurations differ in the binding energy of C20 
fuUerenes in a crystal, AE — E (C^o) — E (Cao)^ where N is the number of fuUerenes 
{N ^ 1 in the macroscopic system). The configuration with the maximum AE value is the most 
energetically favorable one (it has the minimum total energy) and the remaining configurations 
are metastable with respect to this configuration. Moreover, as shown in [13], a crystal of C20 
fuUerenes can have the fee lattice, where C20 fuUerenes are connected with each other by means 
of additional carbon atoms, so that the unit cell consists of 22 atoms. Other three-dimensional 



periodic configurations, which are obtained from various combinations of C20, hkely exist (the 
high chemical activity of C20 fuUerenes can lead to the formation of amorphous structures in 
them) . 

Note that even in the simplest (C2o)n systems, quasi-one-dimensional chains, the number of 
various types of intercluster bonds is quite large [7,10,11]. In these systems, the A£' value is 
maximal in chains with the open- [2-1-2] bonds and is somewhat smaller in chains with "twisted" 
(according to [7]) bonds. Before the search for the three-dimensional (C2o)n structure with 
the maximum AE' value, it is reasonable to study in detail two-dimensional complexes of C20 
fuUerenes. Although the data on the synthesis of cluster films of C„ with n = 20 — 32 exist [14] 
(but have not yet been confirmed), this problem is poorly investigated theoretically [15]. 

Since the strong interaction between clusters in C20 fuUerite can give rise to the transition to 
another atomic configuration (not even crystalline, but amorphous), in addition to the energy 
benefit of a certain structure, its stability is important. As shown in our previous work [16], an 
isolated C20 fuUerene loses its spheroidal shape through the decay into various energetically 
unfavorable configurations, rather than through its transition to the equilibrium "bowl" isomer 
with a minimum total energy. The chain-cluster structure of quasi-one-dimensional (C2o)]\f 
complexes can disappear due both to the decay of one of the C20 fuUerenes and to the coalescence 
of two neighboring C20 into a C40 cluster (the separation of one of the fuUerenes from the chain 
is also possible) [10, 11]. In this case, the height U of the lowest potential barrier around 
the corresponding metastable configuration decreases from « 5 eV in the C20 fuUerene [16] to 
^ 2.5 eV in the (020)2 dimer [17] and further to 2 eV in the (C2o)n chains with N ^ 1, 
where the main stabihty loss channel is the coalescence of fuUerenes [10,11]. In two-dimensional 
systems, both the stabilization of the metastable state due to an increase in the number of the 
nearest neighbors of each fuUerene and the appearance of new stability loss mechanisms can be 
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expected. 

The aim of this work is to analyze the structure and energy characteristics of two-dimensional 
complexes of C20 fuUerenes, as well as their stability. The tight-binding method [18], which 
includes the contribution of all valence electrons (four per atom) to the total potential energy 
E, is used in the calculations. This method does not require such computational resources as 
ab initio approaches and allows the investigation of systems consisting of several hundreds or 
more atoms. For Ceo and C20 clusters, it provides the interatomic distances, binding energy, 
HOMO-LUMO gap, and other characteristics most of which are in good agreement with the 
experimental data and ab initio calculations [16,17,19-21]. For the (C2o)]\f chains, it gives the 
same sequence in the location of isomers with various intercluster bonds in the /S.E value [10,11] 
as more strict density- functional calculations [7]. 

Making use of free boundary conditions, we investigate quasi-two-dimensional (C2o)mxM 
complexes with M — 2 — h and various types of intercluster bonds. By analogy with the 
(C2o)n chains [10,11], it can be expected that the binding energy /S.E is maximal when C20 
fuUerenes are connected with each other by open-[2+2\ bonds in both horizontal and vertical 
directions (configuration A, Fig. la). However, the AE value appears to be much higher (by 
~ 1 eV/C2o) if intercluster bonds are twisted along one of the directions (configuration B, Fig. 
lb). Configuration B is obtained from configuration A by the shift of the horizontal chains 
with respect to each other so that the intercluster bonds in the vertical direction (i.e., bonds 
between fuUerenes of different chains) become twisted. This configuration can be treated also 
as a configuration where vertical chains with twisted bonds are connected by the open-[2+2] 
bonds. The difference between the binding energies AE in the (C2o)]V[xM complexes with 
configurations A and B increases with the size of the system (i.e., with M). The extrapolation 
of the data obtained for M = 2 - 5 to the limit M 00 yields AE = 8.00 eV/C2o and 
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9.54 eV/C2o for configurations A and B, respectively (see Fig. 2). We also calculated A£' for 
configuration B with periodic boundary conditions. The value AE = (9.45 -r 9.81) eV/C2o 
obtained by extrapolating to the limit M — > oo almost does not differ from the value obtained 
with free boundary conditions (a comparatively large error is due to the nonmonotonic saw- 
tooth dependence of AE on M and the small number of points used for extrapolation) . 

To analyze the thermal stability of the (C2o)]V[xM complexes, we used the molecular dynamics 
method. Configuration B with the combination of the open-[2+2] bonds and twisted intercluster 
bonds (see Fig. lb) was taken as initial. At t — 0, all atoms were assigned with random 
velocities and displacements so that the total momentum and total angular momentum were 
zero. Classical Newtons equations of motion were solved numerically with a time step of 
to = 2.72 X 10~^^ s. The total energy of the system (the sum of the potential and kinetic energies) 
remains unchanged in the simulation process, corresponding to a microcanonical ensemble, and 
its dynamic temperature T, which is a measure of the energy of the relative atomic motion, 
was calculated by the formula [22,23] (-Ej^jn) = ^k-^T {3m — 6), where (-Ej^jn) is the time- 
averaged kinetic energy of the entire system, m — 20M^ is the number of atoms in it, and /cg 
is the Boltzmann constant. Since the lifetime of the (C2o)]VIxM complexes to the instant of loss 
of their periodic cluster structure increases exponentially with a decrease in T and the time 
necessary for the numerical calculations at a given T value is proportional to M^, we restricted 
ourselves to the systems with M = 2 and 3 at T = (1800 - 3300) K. 

Sometimes we observed transitions from configuration B to configuration A and vice versa. 
This is due to the nearly equal values of the energies of these configurations (see Fig. 2) and 
likely to a quite low potential barrier between them (according to the results obtained for the 
(C2o)n chains [10], its height is expected to be about 1 eV). However, although such transitions 
are rather fast (occur in a time of 0.1-1 ps corresponding to ~ 1000 molecular dynamics steps), 

5 



they are sufficiently rare and do not lead to the breakdown of the spheroidal shape of C20 
fuUerenes or to the breaking of the periodicity in their arrangement. 

As in the simulation of the (C2o)]\f chains [11], we observed two fundamentally different sta- 
bility loss channels for the(C2o)]\/[xM- decay of one of the C20 fuUerenes and coalescence 
of two neighboring C20 fuUerenes into the C40 cluster. The decay and coalescence are accom- 
panied by a decrease and an increase in the binding energy by (0.5 — 1) eV respectively; 
i.e., the system transits to the energetically less favorable and more favorable configurations, 
respectively. In both cases, the first decay /coalescence act is followed by the decay and/or 
coalescence of the remaining fuUerenes of the complex (see Fig. 3). In particular, it is some- 
times difficult to distinguish whether the sequential coalescence of two C20 fuUerenes into the 
C40 cluster occurs followed by its coalescence with another C20 fuUerene into the Cgo cluster or 
three C20 fuUerenes immediately fuse into the Ceo- cluster. For the same reason, it is sometimes 
impossible to clearly determine which event - decay or coalescence - occurs first initiating the 
integrity loss of the entire system. Nevertheless, the temperature of the (C2o)]V[xM complex 
correlates with the stabihty loss mechanism: at 1800 <T < 2400 K, the periodic cluster struc- 
ture disappears due only to the coalescence of C20 fuUerenes; whereas at T > 2400 K, due both 
to the coalescence of C20 fuUerenes and to their decay, although the decays occur much rarely 
(see Fig. 4). Approximately the same processes occur in the (C2o)n chains [11]; however, in 
contrast to the chains, the separation of one C20 fuUerene from the (C2o)]VIxM complex is not 
observed probably because each C20 fuUerene has a covalent bond with no less than two other 
fuUerenes. 

To determine the activation energy Ea for the coalescence of the C20 fuUerenes in the (020)2x2 
complex, we analyzed the temperature dependence of the lifetime using the Arrhenius formula 
for the microcanonical ensemble with the finite-heat-bath correction [24,25]. We assumed that 
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T-^(r) = Ae-x.^[-Ea/k^T*), where T* ^ T - EJ2C (the microcanonical heat capacity 
is assumed to be C = (3m — 6) /cg, where m = 20M^; we will show below that Ea/2C ~ 
100K< T at T f« 2000 K). As seen in Fig. 4, the calculated dependence of In (r) on 1/T* 
in the first approximation is given by a straight line (because of the rather large spread of 
the "numerical experiment" data, the temperature dependence of the parameter A cannot be 
reliably determined) . Using the slope of this straight line and the point of its intersection with 
the ordinate axis, we determined the activation energy Ea = (2.45 ± 0.2) eV and the frequency 
factor A = (6 ± 1) X lO^'^'s"-^. The Ea value coincides within the errors with Ea = (2.7 ± 0.3) 
eV for the coalescence of C20 fuUerenes into the (020)2 dimer [21]. This value is generally 
averaged over the different coalescence channels. Indeed, analysis of the potential energy of 
the (020)2x2 complex as a function of the coordinates of constituent atoms shows that there 
are many different paths for the transition from configuration B to the configurations with the 
fusedfuUerenes and that the heights U of the potential barriers on these paths vary in a quite 
wide range of C/ = (2.4 - 3.0) eV (we do not exclude the possibihty that barriers with a lower 
U value exist, but we have not found them due to the complex shape of the potential relief). 

In order to collect the simulation results sufficient for determining Ea in the (020)3x3 complex, 
a very long computer time is required. Nevertheless, the comparison of r values obtained for 
certain temperatures in the range T = (2000 — 3000) K with the r (T) dependence for the 
coalescence of the O20 fuUerenes in the (020)2x2 complex (see Fig. 4) shows that the stabihty 
of the system remains almost unchanged with an increase in the system size: Ea 2.5 eV for 
both M — 2 and 3. The complete absence of the decay processes at T < 2400 K indicates that 
their activation energy is quite high (likely, Ea — (3.5 — 4) eV, as for the (020)2 dimer [21]). 

Summarizing the results, we can state that the main channel of the breakdown of the pe- 
riodic cluster structure of the two-dimensional (02o)]V[xM systems is the coalescence of their 
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constituent C20 fuUerenes. The average activation energy of such processes is about 2.5 eV; 
i.e., the hfetime of these systems at room temperature is macroscopicaUy long. Since the sta- 
bihty of the two-dimensional complexes of the C20 fuUerenes is no worse than the stability of 
one-dimensional chains [11], it can be expected that the three-dimensional systems of the C20 
clusters are also characterized by a quite high stability. 

To conclude, we note that the density of the single-particle electronic states at the Fermi 
level in two-dimensional (C2o)|vIxM complexes is finite (see Fig. 5) and depends weakly on 
M for M = 3 — 5. For this reason, it can be expected that, if films of C20 fuUerenes were 
synthesized, they will be conducting. To manufacture such films, one can try to modify the 
method of low-energy beam deposition of neutral clusters [14]. In this case, it is very important 
to select the appropriate substrate material (e.g.. Ceo monolayers were recently grown on Au 
substrates [26]). It is also of interest to perform the numerical simulation of the coalescence of 
initially isolated C20 fuUerenes into two-dimensional complexes. 
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Fig. la. Two-dimensional (C2o)4x4 structures with the open-[2+2] intercluster bonds, config- 
uration A. The binding energy of the fuUerenes is AE — 6.11 eV/C2o- 
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Fig. lb. Two-dimensional (€20)4x4 structures with the combination of the open- [2+2] and 
twisted intercluster bonds, configuration B. The binding energy of the fuUerenes is A£' = 7.08 
eV/C2o. 
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Fig. 2. Binding energiesAE' of the C20 fuUerenes in the two-dimensional (C2o)]V[xM structures 
with various types of intercluster bonds versus 1/M for configurations A (triangles) and B 
(circles). The solid lines are the least-squares approximations. 
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Fig. 3a. Configuration of (C2o)3x3 formed after the decay of one of the C20 fuUerenes. The 
initial temperature is T — 2500 K, the hfetime with respect to the decay is r — 2.7 ps, and the 
binding energy (after relaxation) is AE — 5.93 eV/C2o- 
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Fig. 3b. Configuration of (C2o)3x3 formed after the coalescence of two C20 fuUerenes into the 
C40 cluster. The initial temperature is T = 2220 K, the hfetime with respect to the coalescence 
is T = 59.6 ps, and the binding energy (after relaxation) is A£' = 7.01 eV/C2o- 
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Fig.4. Logarithm of the hfetime (in seconds) of the (C2o)]\/[xM complexes with M — 2 (closed 
symbols) and 3 (open symbols) to the instant of the decay of one of the C20 fuUerenes (squares) 
and the coalescence of two C20 fuUerenes into the C40 cluster (triangles) versus the inverse 
initial temperature T* with account for the finite- heat-bath corrections (see the main text). 
The solid line is the linear approximation for the coalescence processes at M = 2. 
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Fig. 5. Density of the single-particle electronic states for the (C2o)4x4 complex with the 
combination of the open- [2+2] and twisted bonds (see Fig. lb). The energy is measured from 
the Fermi level. 
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